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Abstract. We describe a general algorithm for generating various families of ribbon tableaux and 
computing their spin polynomials. This algorithm is derived from a new matricial coding. An 
advantage of this new notation lies in the fact that it permits one to generate ribbon tableaux with 
skew shapes. 



1. Introduction 

Ribbon tableaux are planar structures generalizing Young tableaux (see [5l [lU] for the classical 
case). These are tilings of Ferrers's diagram by ribbons (diagrams with special shape) labelled with 
integers verifying some vertical and horizontal monotonicity conditions. 

Standard ribbon tableaux (all labels different) have first been introduced by Stanton and White 
in 1985, in order to explain some combinatorial properties of colored permutations Semi- 
standard ribbon tableaux (repeated labels are allowed) go back to the work of Lascoux, Leclerc 
and Thibon [8]. These authors were motivated by the introduction of g-analogues of certain combi- 
natorial identities, and in particular of the famous Littlewood-Richardson rule describing products 
of Schur functions. They obtained g-analogues of decomposition coefficients for any product of 
Schur functions and many questions about these g-coefficients are still open. 

Studying ribbon tableaux is quite a difficult subject, which mainly uses huge numerical ex- 
perimentations. This is why we are interested in finding efficient algorithms for generating and 
computing statistics on them. 

The matricial coding of ribbon tableaux used for numerical experimentations in [HI [9] is deficient 
because it does not give some elementary properties (shape and position of the head of the ribbons 
for example) without additional computations and cannot be generalized to skew shapes. The 
algorithm used at the time was not published and appeared only as a programming example with 
Maple/ ACE in [6] and a distributed version is described in [T3] . 

The aim of this paper is to present a more general algorithm for generating ribbon tableaux and 
computing spin polynomials, using a different and more transparent coding for ribbon tableaux. 
These algorithms are implemented in the combinatorial library MuPAD-Combinat^ which can be 
downloaded at http://mupad-combinat.sourceforge.net/. 



2. Basic definitions on ribbon tableaux 

We will mainly follow [H [5l [TO] for classical notions on partitions and tableaux and [lOj for 
notations related to partitions. Let A and fj, be two partitions such that the diagram of A contains 
the diagram of The skew partition of shape A//u can be defined as the set-theoretic difference 
A — ^ (A is called the outer partition and fi the inner). 
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h(R) 



Figure 1. An 11-ribbon of height h{R) = 4 



Definition 1. Let k be a nonnegative integer. A A; -ribbon R is a connected skew diagram with k 
cells which does not contains a 2x2 square. The first (north-west) cell is called the head and the 
last one (south-east) the tail. The spin is defined as sp(R) = . 



Figure 2. A tihng with 3-ribbons of the skew partition (8, 7, 6, 5, 1, l)/(3, 3, 1) of 
spin 3 

We will denote by Pavk{X/ fi) the set of /c-ribbon tilings of the skew shape A//i. The spin sp{P) 
of a tiling P is the sum of the spins of its ribbons, and the cospin is the associated co-statistic into 
PavkiX/n), i.e: 

cosp{P) = max{sp(J7), U G PaVk{X/iJ.)} - sp{P). 
Ribbon tableaux are labelled tilings verifying monotonicity conditions similar to these of Young 
tableaux. The spin of a fc-ribbon tableau is the spin of the underlying tiling. 



2 3 




Figure 3. A 3-ribbon tableau of shape (8, 7,6, 5, 1, l)/(3, 3, 1), weight (2,2,2,1), 
and spin 3 



Definition 2. A k-ribbon tableau of skew shape X/ ji is a tiling of the skew shape X/ ji by labelled 
k -ribbons such that the head of a ribbon labelled i must not be on the right of a ribbon labelled j > i 
and its tail must not be on the top of a ribbon labelled j > i. The weight of a k-ribbon tableau is 
the vector v such that Vi is the number of k-ribbons labelled i. 

We denote by Tabfe(A//x, i/) the set of all semi-standard /c-ribbon tableaux of shape A//x and 
weight V. 

Definition 3. The spin and cospin polynomials associated to the set Tabk{X/ ji,u) are: 
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If we write sp* = max{sp(T), T € Tabi:{X/ fi)}, the following property holds: 



Example 1. In Tab3{{8, 7, 6, 5, 1), (3, 3, 2, 1)), these two polynomials are: 



sp* Mk) j\ 



G 



(3) 

(8,7,6,5,1), (3,3,2,1) 
(3) 



Aq) = + 17^3 ^ 33^4 ^ 3^^5 ^ ^g^6 ^ 5^7^ 



G 



(8,7,6,5,1), (3,3,2,1) 



i)(g) = + 17/ + 33g3 + 31^2 ^ ^ 5^ 



3. A GENERATING ALGORITHM 

In this section, we describe a new algorithm for generating all the ribbon tableau in Tabfc(A/^, 5) 
and computing their spin polynomials. The main basic idea is to apply recursively the algorithm 
of removing fe-ribbon strips from a partition. 



3.1. A new coding for ribbon tableaux. In this subsection, we extend the coding of [13] which 
is not well adapted to ribbon tableaux with skew shape \/ (j, and does not give immediate access 
to the shape of the tableau. If the partition A is (Ai, . . . , A„), our coding is an n x Ai-array 
defined as follow: 



1- A 

2- A. 

3- A. 



-1 if e ^l, 

p if there is the head of a ribbon labelled p in the cell («, j) 
otherwise. 



By construction, we can immediately read the shape of the corresponding ribbon tableau, and the 
length of the ribbons is obtained by dividing the number of non negative cells by the number of 
positive cells. This coding makes sense when k = 1 because we obtain the classical representation 
of a skew Young tableau (each ribbon is reduced to its head). 

Example 2. The new coding of a 3-rihhon tableau of shape (8, 7, 6, 5, 1, 1) and weight (3, 3, 2, 1).' 

/ 2 \ 

2 3 

2 4 

110 10 

\ooooo3ooy 

The structure used in [13] was a matrix where the entry {i,j) is equal to k if there is the tail of 
a ribbon labelled i in the j'-th column of the shape. We show that the main difficulty in carrying 
out this construction is that we cannot read the shape of the partition without reconstructing the 
entire ribbon tableau. Furthermore, this coding become non injective when we try to generalize 
it on ribbon tableaux of skew shape. With the previous example we would obtain the following 
matrix: 



2 










2 3 








2 
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1 


1 




n 




3 



/03303000\ 
30303000 
00003003 
\00000030y 
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Example 3. The new coding of a 3-ribbon tableau of skew shape (8, 7, 6, 5, 1, 1)/ (3, 3, 1) and weight 
(3,2,2,1); 

/ 1 \ 
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Decoding a fc-ribbon tableau from an array {^ij) is as follow: for each label p considered in 
decreasing order, we fix vertically by the head a ribbon labelled p in the cell {i,j) of the shape if 
Aij = p, and we drive them to the frontier of the partition. 

Algorithm 1 (From coding to ribbon tableaux). 

• Input; an array M = {rriij). 

• Initialize; 

— A/^ < — the shape of the tableau coded by M , 

— u < — the weight of the tableau. 

• For each label i in v: 

— fix ribbons labelled i by the head in each position corresponding to the 
cells labelled i in M , 

— drive these ribbons on the frontier of A, 

— A < — A without previous ribbons labelled i. 

• Output; The k-ribbon tableau corresponding to the array M . 




Figure 4. Running of the algorithm on the previous 3-ribbon tableau of shape 
(8,7,6,5,1,1)7(3,3,1) and weight (2,2,2,1) 
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3.2. Adding and removing fc-ribbon strip. In order to explain our generating algorithm, we 
begin with a general algorithm for adding or removing a fc-ribbon strip from a partition (these two 
operators on partitions come from the representation theory of the quantum algebra Uq{sln), see 
[2j for more details). 

Definition 4. A skew tiling by k-ribbons Q where the tail of each ribbon is not on top of an other 
ribbon is called a /c -ribbon strip. The weight of Q is the number of ribbons in the tiling. Let Q-^ 
(resp. 0jj be the the horizontal strip made of the top cells (resp. the bottom cells) of the columns 
ofQ. 




Figure 5. A horizontal 5-ribbon strip of weight 4 and spin 3/2 

In a ribbon strip, the head of each ribbon lies in and the tails in G|. As a fc-ribbon strip 
has a unique tiling by /c-ribbons, it is completely defined by the positions of all the ribbons's heads 
and the outer partition or all the tails and the inner partition. For adding a fc-ribbon strip to 
a partition, we represent this strip by the integer vector V = (pi, • • • ,Pn) with pi = k \i Q has 
a ribbon with tail in the i-th column of 0|. Similarly, for removing a A;-ribbon strip from a 
partition, the vector V is now defined by pi = —k \i Q has a ribbon with head in the i-ih. column 

of e|. 

Algorithm 2 (Algorithm for adding a A:-ribbon strip to a partition). 

• Input.- the partition A = (Ai, . . . , A„) and the positions vector V = {pi, . . . ,Pr). 

• Initialize.- 

— 6 < — {max{n, r) — 1, . . . , 0), 

— A < — conjugate partition of A. 

• Xi — X + 6 + V. 

• Sorting A: 

- A ^ a{X), 

— I < — the inversions number of a. 

• if X — 6 is a partition: 

- then A ^ conjugate (A - 6) and / ^ {k-i){pr+^..+pr)-l ^ 

- else A < — FAIL. 

• Output.- (A,/) if exists, FAIL otherwise. 

In the previous algorithm, a corresponds to a permutation which permits to sort the vec- 
tor A. Let defined the number of inversions of the permutation a by the cardinality of the set 
{{i,j) such that i < j and a{i) > cr{j)}. 

3.3. Generating algorithm and computation of the spin and cospin polynomials for 
ribbon tableaux of a given shape and weight. We will generalize the generating algorithm 
of ^3j to the case of fc-ribbon tableaux of skew shapes. A basic remark consists in the fact that 
the ribbons labelled i form a k ribbon strip 0j of weight ui. That's why we search recursively all 
the fc-ribbon strip of weight Ui contained in the shape and removable from the outer partition and 
we fill an array with the positions of the head in Ojf. As the algorithm 1 also returns the spin of 
the added/removed fc-ribbon strip, by keeping the spin in each step of the construction we obtain 
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finally the spin of each tableau. In other words this algorithm permits one to compute the spin 
and cospin polynomials without additional computation. 

Algorithm 3 (Generating algorithm for ribbon tableaux and spin polynomials). 

• Input.- The shape X/ fi, the weight v = (z/i, . . . , Up) and k. 

• Initialize: P = and L < — {(T, A,0)} where T is an array filled with -1 in 
cells corresponding to and otherwise. 

• For each weight j from p down to 1, 

For each (T, At, spx) in L and each permutation C of the vector (0, . . . ,0,k, . . . ,k) 

Xi—i'j times 

corresponding to a valid removed /c-ribbon strip, 

— fill in T the cells corresponding to the frontier of At and the non-zero 
coordinates of C, 

\t < — At without the A;-ribbon strip corresponding to C, 

— spT * — spT+ the spin of the previous fc-ribbon strip, 

— if C At then L < — L U (T, At, spr)- 

• For each tableau T in L, P< — P + q^P^ . 

• Output.- 

— The list L of all the k-ribbon tableaux of skew shape \/ fi and weight v, 

— P the spin polynomial of this set. 

By studying the progress of this algorithm with huge set of ribbon tableaux we remark that at the 

bottom of the tree there are a lot of nodes, but in there is only few different shape for the remaining 
partition. For example the number of nodes for A = (6, 6, 6, 6, 6, 6) and v = (3, 1, 1, 1, 1, 1, 1, 1, 1, 1) 
is the sequence 3, 12, 48, 198, 780, 2940, 10080, 31080, 81480 and finally 43680 but the number 
of different remaining shape are 3. 9, 16, 27, 33, 38, 33, 27, 16, 1. Our previous algorithm search, 
at step i, all the possibility to retire a k ribbon strip from the remaining partition, but in fact 
we need to test all the possibility only on the few remaining partitions. That's why a recursive 
implementation of the previous algorithm with a remember option seems to be the most efficient 
way to generate ribbon tableaux. 

Example 4. In the case of 3-rihhons with A = (9, 9, 9, 9, 9, 9, 9, 9, 9) and standard weight jj, = (1^^), 
we have the following spin polynomial: 

414315330 + 8286306600^ + 85027356570(^2 + 588666753870^'^ + 3062543589300g^ + 
12659483135520^5 + 42941179272810(?6 + 1219126827839709^^ + 293410572110760^^ + 
603798294330270^9 + 1068859924958280gi° + 1634693172838050gi^ + 
2166452577489720^12 + 2492870571244950^^3 + 2492870571244950^^^ + 
2166452577489720gi5 + 1634693172838050^1^ + 1068859924958280gi'^ + 
603798294330270gi^ +293410572110760^19 + 121912682783970^20 _^4294ii79272810g2i + 
12659483135520g22 + 3062543589300^^3 + 588666753870^24 _^ 85027356570g25 + 
8286306600^26 + 414315330^27 

which correspond to a computation over 16 882 686 792 972 000 ribbon tableaux. 
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